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Abstract 



We give a procedure for reconstructing a magnetic held and elec- 
tric potential from boundary measurements given by the Dirichlet to 
Neumann map for the magnetic Schrodinger operator in R™, n > 3. 
The magnetic potential is assumed to be continuous with L°° diver- 
gence and zero boundary values. The method is based on semiclassical 
pseudodifferential calculus and the construction of complex geometri- 
cal optics solutions in weighted Sobolev spaces. 

1 Introduction 

Let Q C R n , n > 3, be a bounded domain with C > boundary. We consider 
the magnetic Schrodinger operator 



where Dj = i^f-, W G L°°(0;C") is the magnetic potential, and q G 
L°°(Q; C) is the electric potential (the coefficients can be complex valued). 
Assuming is not a Dirichlet eigenvalue of i?w,g i n ^> the problem 



it 



Hw. 




in Q, 
on dQ 



has a unique solution u = Uf G for any / G H 1 ^ 2 (dfl). 



The boundary measurements are given by the Dirichlet to Neumann map 
(DN map), defined formally by 



(A w>q f,g) = / {Vu f -Ve g + iW ■ (u f Ve g - e g Vu f ) + (W 2 + q)u f e g )dx 



where Vh E H 1 ^) solves the adjoint problem H_w,q v h = in f2 with v^qq = 
h, and is any H 1 ^) function with eh\dn = h. Then A\y, q is a bounded 



The gauge transformation W i— > W + Vp, where p E V^ 1,oo (0; C), trans- 
forms the magnetic potential to a gauge equivalent potential but preserves 
the magnetic field curlVF. If additionally p\qq = then Aw/ + vp, g = A\y, q , 
which means that boundary measurements are preserved in gauge transfor- 
mations which respect the boundary. 

We are interested in recovering curl TV and q from Aw >q . This is a typ- 
ical inverse problem where one wishes to know the interior properties of a 
medium by making measurements at the boundary. It is related to the ex- 
tensively studied inverse conductivity problem of Calderon [3]. In fact, most 
known results in dimensions n > 3 (starting in Sylvester-Uhlmann for 
later work see (SO]) reduce that problem to recovering q from Ao )9 . 

The problem is closely related to inverse scattering at fixed energy, which 
was studied for the magnetic case in Eskin- Ralston See also Sun [21] 
for an application of the method in One motivation for the present 
study has been to understand the approach of [S| and to clarify the role of 
pseudodifferential operators in inverse problems for first order perturbations 
of the Laplacian. 

Previous results for this inverse problem concern unique determination 
of the coefficients, and they state that Ayp 1)9l = Ajy 2j(?2 implies curlTVi = 
curl W% and q\ = q2 in $7, under varying assumptions on Wj, qj, and 0. Sun 
|23| proved this in the case where the curlPFj are small. Panchenko |19j 
proves a similar result in a less regular setting. The smahness assumption 
was removed by Nakamura-Sun-Uhlmann ^Sj who considered C°° coeffi- 
cients. Based on the method in |15| . the smoothness assumption was reduced 




Here v is the outer unit normal to d^l. More precisely, if /, g E H 1 / 2 (9f2) 
we define Ay^ >q using the equivalent weak formulations 




Vv g + iW ■ (e f Vv g - v g Ve f ) + (W 2 + q)e f v g ) dx 



map H x l 2 (dQ) -»■ H'^dSl). 
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to C 1 by Tolmasky |27| , and to Dini continuous in [2Uj . All these results are 
nonconstructive. We also mention the recent papers |2J, jS] and which 
consider boundary determination, partial Cauchy data and stability for this 
inverse problem. 

In this paper we give a constructive algorithm for recovering curl W and 
q from A\y^ q . Such an algorithm is well known in the case W = and is 
due to Nachman |13| and also Novikov ^H] . This result combines ideas from 
scattering theory with so called complex geometrical optics (CGO) solutions 
of the Schrodinger equation, which were introduced in the fundamental pa- 
per [23]. The main point is that the CGO solutions are defined globally 
and are unique in a certain sense. The problem in extending the results 
of Nachman to nonzero W has been that the main method of producing 
CGO solutions in that case, the pseudodifferential conjugation technique of 
Nakamura-Uhlmann [IJj], only works in bounded domains and there does 
not seem to be a proper notion of uniqueness for solutions. 

We will give a global version of the Nakamura-Uhlmann technique which 
will produce global CGO solutions with uniqueness in the proper weighted 
Sobolev spaces. To do this we apply semiclassical pseudodifferential cal- 
culus. This is largely equivalent with the parameter-dependent calculus 
used earlier in such results, but it simplifies the proofs. A main new ele- 
ment in our approach is a variant of the pseudodifferential cutoff technique 
used by Takeuchi |26| and Kenig-Ponce-Vega jHj in the context of nonlinear 
Schrodinger equations. We remark that semiclassical notation was also re- 
cently used by Kenig-Sjostrand-Uhlmann ^U] who studied the problem of 
recovering coefficients from partial boundary measurements. 

We record some notation. Let A = J^-D 2 , A^ = A + 2( ■ D, and 
D c = D + C, where D = (Di, ... , D n ) and ( e C n , ( 2 :=(■( = 0. The 
^-dependent operators A^ and arise naturally in the construction of 
CGO solutions. For 5 E R we also use weighted I? spaces L 2 with norm 
II/IIl 2 = II ( x ) <5 /IIl 2 where (x) = (1 + Ix) 2 ) 1 / 2 , and weighted Sobolev spaces 

5 

H$ with norm ||/||_f/| = IK^^/ll-f/ 8 - If ^ is a function space we write X c for 
the set of compactly supported functions in X, and Xq for all functions in 
X supported in Q. 

The construction of CGO solutions follows from the norm estimates in 
the following theorem. 

Theorem 1.1. Let W E C c (R n ;C n ), q E L^°(R n ;C), and -1 < 5 < 0. If 
( E C n with £ 2 = and \C\ is large enough, then for any / E L 2 , 1 (R n ) the 
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equation 

(A c + 2W ■ D c + q)u = f (1) 

has a unique solution u G ff](R n ). Furthermore, u G -ff|(R ra ), and u 
satisfies for < s < 2 

\\u\\ H ; ^ClCrHfh^, 
where C is independent of £ and /. 

Theorem ll.ll generalizes results in |25j to operators with first order terms. 
Using the global CGO solutions to Hy/, q u = obtained from this theorem, 
we may extend the results of Nachman to obtain a constructive algorithm 
for recovering curlVT and q from Ky/ q . The main result is as follows. 

Theorem 1.2. Let C R n , n > 3, be a bounded simply connected C > 
domain. Suppose W G Cn(R n ; C n ) with D ■ W G L°°(R n ; C), and suppose 
<7 G L^(R n ; C). Also suppose that is not a Dirichlet eigenvalue of Hyy jq 
in £1. Then determines curl TV uniquely and constructively. Further, if 
W is C 1+e and is C 2+e for some e > 0, then one may construct q from 

The reconstruction procedure for curl W is outlined in the following four 
steps. 

1. From the knowledge of Ayy tq , one may determine the boundary values 
u (\dn °f a CGO solution as the unique solution of a boundary 
integral equation on dQ. 

2. From A\y jQ and u^qq one computes a scattering transform tw,q(£,,0- 

3. The expression Rw, q (£,, m) = lim^oo s _1 tvF,q(?) s/i) is essentially the 
Fourier transform of curll^. 

4. curlVT may be computed from Rw,q using the inverse Fourier trans- 
form. 

The structure of the paper is as follows. Section 2 contains some facts 
on semiclassical pseudodifferential calculus, and Section 3 contains estimates 
for d equations which will be needed later. In Section 4 we prove Theorem 
11.11 where the main step is to conjugate the first order term into a lower 
order one using pseudodifferential operators. Section 5 discusses equivalent 
problems which characterize the CGO solutions. In Sections 6 and 7 we 
reconstruct the magnetic field and electric potential, respectively. 
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2 Semiclassical pseudodifferential calculus 

Our method will involve pseudodifferential operators depending on a small 
parameter h. We will collect the required properties of these operators here. 
See [I] for more details. 

Definition. If < a < 1/2 and m G R, we let S™(R n ) be the space of 
all functions a(x,£;h) where x,£ G R n and h G (0,/io], ho < 1, such that 
a(-;h) G C°°(R 2n ) and 

\d^a(x,^;h)\ <C aP hr^\{i) m 
for all a, j3. If a £ S 1 ™ we define an operator A = Op h (a) = a(x, hD) by 

Af(x) = (2tt)-" / e ix <a(x, h£; h)fe) d£. 

The class of such operators is denoted by OpS™. 

Note that we use the standard quantization instead of Weyl quantization 
in the definition of the operators. 

Proposition 2.1. @] Let a G S™ with m G R and < a < 1/2. 

(a) Op fe (a) is a continuous map 5? — » and oS^' — > =5^'. 

(b) If m = then Op^(a) is bounded L 2 — > L 2 , and there is a constant C 
with 

||Op h (o)|| L2 ^ < C 

for < h < ho. 

(c) d Xj Op h (a) = Op h ia)d Xj +Op h (j% 

(d) The adjoint Op h (a)* = Op h (a*), where a* G S™ satisfies for any N 

h\ a \d?D a a 
a*=Y, - + h N ^S?. 

\a\<N 

(e) If a G S™ and 6 G Sf then Op h (a)Op h (6) = Op h (c) where c G S™ W 
satisfies for any N 

c= ^ + h N(l-2a) gm+m 

\a\<N a ' 
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Also, [Op ft (o),Op h (6)] = Op h (d) where d £ S™ +m ' and 

d = -H a b + h 2 ^- 2 ^S^ +m ' 
i 

where H a = V^a • V x — V x a ■ Vg is the Hamilton vector field of a. 

We will also need the boundedness of Op S® operators on weighted 
Sobolev spaces. Here it is natural to use semiclassical spaces Hg h , defined 
by ||/||fff, = \\{hD) s (x) 5 f\\ L 2 . If s > is an integer an equivalent norm is 

given by E H <^ H P q /IIl1- 

Proposition 2.2. Let a £ S® with < a < 1/2 and s ,5o > 0. Then 
Op h (a) is bounded H% h (R n ) — > Hg h (R n ) for any s,(J £ R, and there is a 
constant C with ||Op^(a)||ij-| h -+Ht h — C whenever \s\ < sq, \5\ < So, and 
< h < h . 

Proof. All constants below are independent of h. We begin by showing that 

\\(x) s (hDyAf\\ L 2 <C\\(x) s (hDyf\\ L2 . (2) 

It is enough to take s = since otherwise we may consider {hD) s A{hD)~ s , 
which is in Op S® by Proposition If Tn is a nonnegative integer we define 
Tf(x) = (x)~ 2m A((x) 2m f). For / G y one has 

Tf(x) = (27T)-™ J e™<a{x, K){x)~ 2m (I - ^TfiO 

= (2tt)-" J (I- A 5 r«x)- 2 "V^a(x, hZ))fe) d£. 

It follows by differentiation that T is in OpS®, and Proposition 12.11 gives 
(|2"jl for 5 = —2m. The estimate for <5 < follows from the Stein- Weiss 
interpolation theorem, and for <5 > by duality using the fact that Op h (a)* 
is a pseudodifferential operator. 

It remains to relate (J2| to the norm \\(hD) s (x) s f\\x J 2. If k > is an 
integer then || (hD) 2k (x) 5 f\\ L 2 < C£ |Q| < 2 J(x} a (W})<7|| L2 , using h < 1. 
We claim that 

IK^^^r/IU^cii^^^/ll^, H<2fc. (3) 

In fact, ||(x) 2m (/iD) a /|| L 2 = C||(L> ? ) 2m r(e)(/i0 2fc /1lL2 ifm > is an integer, 
where r(£) has bounded derivatives of all orders. Differentiation gives that 
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this is bounded by C]C|/3|<2m 

\\D^{hi) 2k f\\ L 2. Going back to the x-side 
gives (J3J) for 5 = 2m, and the estimate follows for 5 G R by interpolation 
and duality. Using (J3J) implies that 

\\(hD) 2k (x) s f\\ L2 <C\\(x) s (hD) 2k f\\ L2 . 



The last estimate applied to Af and then (J2J) give 



UhD^W'AfWvKCUxyihD^fh 



|a|<2fc 

||(W?) a <*>'/ll^ 



|a|<2fc 

<C\\(hD) 2k (x) 5 f\\ 



L 2 ■ 



The interpolation (now the Stein- Weiss theorem on the Fourier side) and 
duality give the desired result. □ 

3 Estimates for d equations 

In this section we collect some elementary estimates for equations of d type 
in R n . Let /j, = 71 + 272 where jj G R n , |7j| = 1, and 71 • 72 = 0. The 
operator = fj, ■ V is just d Xl + id X2 in different coordinates, so it has an 
inverse given by 

If 1 

N Z X fi x ) = 7T ~—f( x - 2/i7i - 2/272) dyi dy 2 - 

2vr y R2 yi + iy 2 

This operator satisfies the following. 

Lemma 3.1. Let / G TU fc '°°(R n ) with / = for \x\ > M. Then u = 
N^ 1 f G PU fe,00 (R n ) solves the equation N^u = f in R n and satisfies for 
\a\ < 

|«9 Q u(x)| <C(M)||a Q /|| LO c(x T )- 1 Xfl(o,Af)(a;±) (4) 

where is the projection of x to the plane T = span{7i, 72}, xj_ = x — xt, 
and Xs(o,Af) * s the characteristic function of -6(0, M). 

Proof. Since the statements are rotation invariant we may assume 7,- = ej 
(the jth coordinate vector) and iV^ = d = d Xl + ic?^. Then u = iV" 1 f is 

<x) = ±- f —^—f(x' - y', x") dy>. (5) 
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We have u G L°° since / € L£°. If y> € C°° satisfies = for \x\ > R we 
have 

(u,-dcp) = -— / — — f(x' -y',x")dip(x)dy' dx. 

2?r Jr™ Jr 2 yi + ^2 

The integrand is nonzero only if \y'\ < M + R, which justifies using Fubini's 
theorem. A change of variables and another use of Fubini's theorem gives 

(u, -d V ) = I /(*) ( - -L / -4— ^(x' + y', s") dy') dx. 

The inner integral is <£>(x), which shows that du = f even when / £ L£°. 

It is enough to prove (jlj) for a = 0. From (JSJ) we see that it(x) = for 
| a;" | > M. If | a;' | > 2M then |y'| > |a?'|/2 on the support of the integrand in 

and Q follows. Since Q is easy when < 2M we obtain the required 
result. □ 

We will need a version of Lemma 13. II where / and fx depend on a para- 
meter. Let V Q R n be an open set and let 7j(C) (j = 1, 2) be C°° functions 
of £ G V which satisfy 

l-s< | 7i (C)| <l + e, |ti(0 • 72(01 <e (6) 
and also |<9 Q 7j(£)l < Mi for |a| > 1. 

Lemma 3.2. Let e > be small enough and let f{x,£) G C°°(R n x V) 
satisfy /(x,C) = for |x| > M. Then the function 

If 1 

u ( x ,0 = tt / -—f{x-yili{Q-y2l2{C),i)dyidy2 

2vr J R 2 yi + iy 2 

is in C°°(R n x V), solves (71(C) + 172(C)) " ^ n = / in R n , and satisfies 

\6%d%u(x,g)\ <C Q/3A fM 1 ( J] ll<929|/||Loo (R n x y ) )(xr) l/3hl XiJ(0,M)(2;±y 

h+S\<\a+/3\ 

where xt is the projection of x to the plane T = span{7i(C), 72(C)} an d 
x± = x — XT- 
Proof. Since / is smooth and compactly supported in x it is easy to see that 
u is smooth and solves the given equation. We have 

<%•(/(£ - J/i7i (C) - J/272 (£),£)) = V x /(x - yi7i(C) - 2/272(6,0- 

(-yi%7i(C) - 2/2^.72(0) + d^f(x - yi 71(C) - J/272 (£),£)• 
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Induction and the estimate on the derivatives of 7j(£) imply that 

\P\ , 

\d%dj! u (x,o\ < c aPMl Y, E / \y'\ r ~ l \d2dlf{x- yill {0-y2i2{0^)\dy'. 

|7+<5|<|a+/3| r=0 ^ R2 

(7) 

The integrals in (JJJ) are over the set 

K(x, = {y'eR 2 ;x- yi 7 i(0 - y 2 7 2 (£) G 5(0, M)}. 

We first note that if \x±\ > M then the right hand side of ((7J) is zero, so we 
may assume \x±\ < M. For the behaviour in xt we note that 

y r eK(x,0 e> x-A(y 1 ,y 2 ,0) t £B(0,M) 

(j/i, 2/2, 0)* e A _1 x + yT^O, M), (8) 

where A = .A(£) = (71 (0> • • • >7n(C)) i s a matrix written in terms of col- 
umn vectors, and where 77(C) (3 < j < n) are any orthonormal basis of 

{7i(£), 72(0}^ 

We need to estimate the matrix norm \\A 1 || = sup| x | =1 |^4 1 x\. From 
we obtain ||A|| < 2 for small e, so > |x|/2. An easy calculation 

using the orthogonality properties of the 7,- gives in terms of row vectors 

A' 1 = (ajl + b 7 i cy\ + dji yl 7*) 

where a, b, c, d are obtained from 

a b \ __ 1 / I72I 2 -71 • 72 



c dj I71H72I 2 - (71 -72) 2 V -7i • 72 l7i| 2 

It follows from © that for small e one has a, d ~ 1 and b, c ~ 0, implying 
that || A- 1 1| < 3/2 for small e. 

Suppose I xt I > 12M. If y' € if(a7, then (JHJ) gives 

|y'| > IA^xtI - \A" 1 x±\ - P _1 ||M > \x T \/2 - 3M/2 - 3M/2 > |x T |/4 

and also 

\y'\ < \A~ l x T \ + l^" 1 ^^! + p _1 ||M < 3|x T |/2 + 3M < 2|x T |. 

Since (jHJ) implies that K(x, £) is contained in a ball of radius 3M/2, we obtain 
the desired estimate for u from in the case \xt\ > 12M. If |xt| < 12M 
then (jHJ) implies \y'\ < 21M, and the desired estimate follows also in this 
case. □ 
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4 Proof of Theorem 11.11 

We want to prove the existence, uniqueness, and norm estimates of solutions 
of 

(A c + 2W ■ D c + q)u = f 

with / G L^ +1 (R n ). This will be based on the following fundamental esti- 
mates for the inverse of A^, which imply Theorem ll.ll in the case W = q = 0. 

Proposition 4.1. Let -1 < 5 < and let ( G C n , ( ■ Q = 0, |C| > 1. Then 
for any / G L^ +1 (R n ) the equation 

A c u = f 

has a unique solution u G L^(R n ). The solution operator, denoted by A^" , 
is a bounded map from L^ +l to H$ and satisfies for < s < 2 

IIA^/llHf ^Coicr 1 !!/!!^, 

where Co = Co(n, 8). 



Proof. The main estimate is the case s = 0, which is proved in Since 

1 



we could not find a reference for the Hg result we will give the proof here 



following the s = 1 case in pQ. 

Take 0(f) G C£°(R n ) with = 1 for |f| < 4|C|, = for |f| > 8|C|, and 
|V0j < C7/|C|- For / G L^ +1 we write 

D J D k A- 1 f = T(A- l f) + Sf (9) 

where 

tu = p-HtiSkftQW)}, 



s/ ^- 1{ k^l) /({)} . 

We claim that when — 1 < 8 < 1, 

[|r«[| £ a <C|C| 2 ||n|| L 2. (10) 

For (5 = this follows by Fourier multiplier properties since |£j£fc0(£)| < 
C| CI 2 - For 5 = 1 the statement is equivalent with 

\\Tu\\ Hl <C\(\ 2 \\u\\ H i. 



10 



One has 

v(fu(0) = v(^ fc 0(O)fi(O + ^fe0(OVfi(O- 

Since |V(&&0(£))| < C|CI, we have 
using | C| > 1. This gives (fTU)l for 5 = 
and an interpolation gives for 
For S we claim that 

[IS/Hjj < C||/||^ (11) 

for -1 < 6 < 1. We may write 5/ = ^~ x {i/)f} where = ^fec-f }) ' If 
|£| > 4|C| then | |e| 2 + 2C - ^1 > |£| 2 - 2|C| \Z\ > |£| 2 /2, so |^| < C and | V^| < 
C/|£|. Similar computations as for T imply pip for —1 < 5 < 1. Then 0, 
the s = estimate for A^ 1 , (fTU|) . pip, and the embedding Lg +1 —> Lg imply 

||^D fc A^/ll^<C|C|[|/||^ +1 

for — 1 < 8 < and |£| > 1. This gives the s = 2 estimate for A^" , and an 
interpolation gives the estimate for < s < 2. □ 

An easy perturbation argument using Proposition 14.11 proves Theorem 
ITTl in the case where W = 0. The perturbation argument fails in the case 
where W is nonzero and large. Following Nakamura and Uhlmann \^g^ we 
will use pseudo differential operators to conjugate the first order term into a 
zero order term, so that the perturbation argument can be applied. 

First we will write the equation in semiclassical notation. If £ € C" 
satisfies £ • £ = 0, then we have ( = fi/h with \i = 71 + 272 where 71 , 72 £ R n 
satisfy |tj| = 1 and 71-72 = 0, and h = y/2/\Q is a small parameter. 
Note that Proposition 14. II gives IIAr 1 /!! H i < Ch\\f\\ L 2 , and one also has 

> S,h (5+1 

\\D c f\\ L 2 < Ch^WfWrri . With this notation, 

s S 5,h 

A c = h- 2 q(hD), 
W -D c = h~ l r(x,hD) 

with </(£) = £ 2 + 2/i • £ and r(x, £) = M^(x) • (£ + /i). To take care of the 
nonsmooth symbol r we will make a /i-dependent decomposition r = r" + r", 
where is a smooth approximation and r 9 is a remainder which will have 
small norm on suitable spaces when h is small. 



\\Tu\\ m <C(|C| 2 +|CI)IH|^ <C|C| 2 ||n||^ 
= 1, duality gives the estimate for 5 = — 1, 
— 1 < 5 < 1. 
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Next, we will show that one may conjugate the first order term r* into a 
zero order term. In the proof we need some facts about the symbol q. Since 
l(0 = (£ + 7i) 2 — 1 + 2«72 • £, this symbol has zero set 

g - 1 (0) = {$€R n ; |^ + 7l | = l,e- 72 = 0}. 

For £>0we will consider the neighborhood 

[7(e) = {£ G R" ; 1 - e < |e + 71 1 < 1 + £, |£ ■ 72 1 < 4- 

Lemma 4.1. Let 0<cro<(J<l/2 and let r" G 5^ (R n ) have the special 
form r*(x,£) = t0(x) • (f + /x), where W« G C°°(R n ; C n ) satisfies 

\d a W*(x)\ < C a h- a °W, 
W*{x) = for |x| > M. 

Then there exist a, b, tq G 5°(R n ) and e = e(o" , a) > such that 

(Q + 2/i J R a )A = 5Q + /i 1+e i? . (12) 

Further, (x)ro G 5° so that i?o is bounded from Lg to L| +1 , and for any 
so, So > there is ho < 1 such that whenever h < ho, A and B are bounded 
and have bounded inverses on Hg h for |s| < so,\5\ < a h with norms 
bounded uniformly in h. 

Proof. If a G S® then a direct computation using the special forms of Q and 
1$ implies that 

(Q + 2hF$)A = AQ + hOp h (^H q a + 2r s a) + h 2 Op h (A x a + 2W 9 ■ D x a) . (13) 

The last term is in h 2 ~ 2a Op h 5°, and looking at (JEJ) we would like the 
middle term to vanish. If a = with <j> G S®, this would mean that 

(£ + m) ■ = -r« 

since = 2(£ + /x) • V x . The operator (£ + 71 + 372) • looks like d Xl + «9x 2 
in different coordinates provided that £ G g -1 (0), but degenerates away 
from g _1 (0). Therefore we will only work in a neighborhood of g _1 (0) and 
introduce a cutoff G C£ (R fl ) with ^ = 1 in U(e/4) and V" = outside 
of U(e/2), with e as in Lemma 13.21 This will give a symbol w(x,£) with the 
following properties. 
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Lemma 4.2. The function 



1 



1 



-^{£)r\x - yi(£ + 71) - y2 r Y2,€)dy 1 dy 2 (14) 



27T J R 2 yi + iy 2 
is C°°, solves the equation 

(e + /x)-V s «; = -^(Or l, (x,0, 
and satisfies the estimates 



(15) 



Proof. This follows from Lemma 13.21 with V = U(e), 7i(£) = ^ + 7i, and 

72(0=72- □ 

Note that w is compactly supported in £ but does not have good be- 
haviour in x. To take care of this we will need another cutoff x( x ) £ C£°(R n ) 
with x = 1 on -B(0, M), and we will define 

4>(x,0 = X(.h e x)w{x,£) 
where 9 = a — o~q. Then Q15JI shows that <f> satisfies 

\d«d^{x : z)\<c a p{x)- l h-°\ a+ P\. 

Thus <p and (x}4> are in 5°, and also a = e 1 ^ is in S®. We have 
-H q a + 2r tt a = 2(f + n) ■ (V^)e^ + 2rh i(f> 



2e 
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Since h < 1 we have r" = x^ 9 ^) 7 "" an d 
1 



-# a + 2r B a = 2e 



icf> 



(1 - ^(0)r* + h e (£ + n) ■ wV X (h e x) 



The second term is compactly supported in x and £, but the first term is of 
first order and we are not yet in the situation of (|12|) . Here we are saved 
by the fact that the main operator Q in Q12j) is elliptic on the support of 
1 — ip(C), and we may define b G S® by 

?(0 
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Then (fT3]l becomes 

(Q + 2hR*)A = BQ + h 1+e O Ph (2e^[(ti + ^ ■ wV X (h 9 x)}) 

+ h 2 - 2a Op h (h 2a [A x a + 2W* ■ D x a]). 

Choosing e = min{#, 1 — 2a} gives (TT21) with r G 5°, and one even has 
(x)ro £ S 1 ^. Proposition 12,21 gives that R$ is bounded from L 2 to with 
norm bounded uniformly in h. 

It remains to show that A and B are bounded and invertible on Hg . for 
|s| < so , 1 5| < when /i is small enough, and that all norms are bounded 
uniformly in h. The boundedness of A and B follows from Proposition 12.21 
To show invertibility we note that 1/a = e~ l< ^ is in S® and 

Op h (a)Op h (l/a) = I + fc^OpfcM 

where m G by Proposition l2.ll Since Op^(m) has bounded norm on Hg h 
for \s\ < so, \5\ < 5o, the operator / + h l ~ 2r7 0p h (m) is invertible on these 
spaces if h is small enough. Then also A is invertible with norm of the inverse 
uniformly bounded in h. The same applies to B since b = a + hOp h Sj . □ 

Remark. Lemma f4.1l is a global nonsmooth version of the pseudodifferential 
conjugation technique in (see also ^Zj). Similar ideas have been used 
in inverse scattering [Sj, jS], nonlinear Schrodinger equations [H] and 
periodic Schrodinger operators [22] . 

The problem in extending the method to the global case is seen in (|15|) . 
where the derivatives in £ of the symbol grow in x. This behaviour leads 
to poor global properties. A solution, presented in and jHj, was to 
multiply a symbol a satisfying \d"d^a(x,^)\ < C a p (x) 1^1 by a cutoff 

x(Ro(x) / {£,})■ The new symbol is of type (0,0), hence bounded on L 2 , and 
the error term which appears in the equation because of this modification 
is of lower order. One can even invert related operators on L 2 by adjusting 
the parameter Rq. 

In the present case there is an additional parameter h which may be taken 
arbitrarily small, and additionally |£| < Ch~ l in the support of w(x,h^). 
Therefore our cutoff has the simpler form x{h e x), where is chosen so that 
the new symbol falls into S® with a < 1/2. 

We proceed to prove the main norm estimates. 

Proof, (of Theorem II .1|) The proof is given in three steps. 
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Step 1 : A decomposition 

Let if 6 C c °°(R n ), ip > 0, if = 1 for \x\ < 1/2, and ip = for |x| > 
1. Let also f> £ {x) = e~ n (p{x/e) be the usual mollifier. We will use the 
decomposition (as in [2U]) 

W = + W b (16) 

where = W * ip e is a smooth approximation of W, and we make the 
specific choice 

e = h a ° 

where < ctq < 1/2. Then W is a small remainder term, and one has 

\d a W*(x)\ < C a h~^\ a \, 
||W b || L oo as /i -» 0, 

the second estimate by the continuity of W. 
Step 2: Existence 

Using the decomposition (fl6|) . we write the equation Q as 

(A c + 2W^ ■ D c + 2PU b • D c + 9 )u = /. (17) 

Choosing a with <jq < a < 1/2 and so = 5o = 2, Lemma 14.11 gives /io < 1 
and a, b, tq € 5° with 

(A ( + 2W^ ■ D C )A = BA ( + h- 1+£ R . (18) 

We will assume h < ho, so A and B will be invertible. We look for a 
solution of (|17|) of the form it = A^ 1 ?; for v G ^i+v Then u = AA~ 1 A^ 1 v, 
and inserting this in (|17[) and using (|18[) gives 

(SA C i4 _1 A^ 1 + /T^^Ro^A" 1 + 2W b ■ D^A' 1 + gA^ 1 )?; = /. 

We will show that in the operator on the left, the last three terms are small 
perturbations of the first term when h is small. 
Consider the operator 

M = BA C A _1 A^ 1 . 

From (|18|) we get 

M = 1 + 2W* ■ D C A- 1 - h' 1+E R A' 1 A7 1 , 
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and therefore M is bounded on Lg , 1 with norm bounded uniformly in h. It 
is easy to see that M has the inverse 

Nr 1 = A^AA^B' 1 . 

Similarly from 1)18|) we obtain 

M~ l = I- 2W S ■ D c AA- l B~ l + h~ 1+£ RoA- 1 B- 1 

which is again bounded on L 2 +1 , with norm bounded uniformly in h. Also, 
using the mapping properties of the related operators and the decay of 
, we have 

Wh-^RoA^A- 1 + 2W h ■ D.A- 1 + qA- l \\j2 _ >T 2 = o(l) 

II (_ > ~5 uly S+l ^S+l v ' 

as h — > 0. Then we obtain a solution u of Q in the form 

u = A^ 1 M- 1 (I+h- 1+£ R A- 1 A^ 1 M- 1 +2W b -D (: A^ 1 M- 1 +qA^ 1 M- 1 )- 1 f. 

Thus u = A^ v with < CII/IIl^ . The norm estimates for u follow 

from Proposition 14.11 

Step 3: Uniqueness 

It is enough to show that if u £ Hg satisfies 

(A c + 2W* ■ D c + 2W b ■ D ( + q)u = 0, (19) 

then u = 0. We use Lemma 14. II It follows that u = Av for v = A~ 1 u £ Hg, 
so that v satisfies 

(BA ( + h- 1+£ R + 2W b ■ D ( A + qA)v = 0. 

Applying B~ l from the left we get 

A ( v = -{h- 1+e B- x RQ + 2B~ x W b ■ D ( A + B^qA^. (20) 

The right hand side of (|2*U)) is in L 2 g +1 since Rq and qA map L 2 S to L 2 S+1 
and W b ■ D^A maps Hg to L 2 +1 . We are now in the situation of Proposition 
14.11 and using the Hg estimate of that Proposition to (|2Uj) implies 

\\v\\ Hj < CQ\(hr 1+e B- 1 R + B-\A)v\\q +i + \\2B~ X W* ■ D Q Av\\ L , +i ). 
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The first term on the right is bounded by C|C| IMI.L 2 and the second term 

S 

is bounded by C[|(a?)W ||£°a (||u||#i + KIIML 2 )) with C independent of £• 

s s 

Choosing \(\ so large that || (aj)!^^ ||z,°° < 1/(2C), we have that the coefficient 

of |H|^i on the right is < 1/2, and we may move this term to the left. We 

s 

are left with the estimate 

\Mh} <C|CINz3 (21) 

with C independent of (. 

Finally, we use the L? estimate of Proposition 14.11 to (|2U|) . This implies 




The first term inside the parentheses is < C | C 1 1 e |Ml£ 2 > an d the second is 
< C\\(x)W b \\ L ™ (\\v\\„i + |C| IMIl 2 )- Using CD this becomes 

5 5 
\\v\\ L 2<C(\C\- £ +\\(x)W^)\\v\\ L 2 

where C is independent of Q. Choosing \Q\ large enough we obtain (HI £2 < 

s 

^\\v\\ L 2, implying v = and also u = 0. □ 



5 Equivalent problems 

In this section we show that complex geometrical optics solutions for the 
magnetic Schrddinger equation can be characterized in several different ways. 
The treatment is almost completely analogous with [T^. We begin by stat- 
ing the main result and explain the notation later as we go along. 

Proposition 5.1. Let f2 C R n be a bounded domain with C 1,1 boundary. 
Suppose W E L~(R n ; C n ) with D ■ W £ L°°, and suppose q 6 L™(R n ; C). 
Also suppose that is not a Dirichlet eigenvalue of Hw q in Q. Let £ G C n 
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with £ 2 = 0, and consider the following four problems: 

H W)q u = in R n 
u = e iC -- x (l + u) with to G A^L%, 

u + G ( * (2W • Du + {W 2 + D ■ W + q)u) = e^ x in R n 

i) An = in W 

ii) u G H 2 (Q' R ) for any R > R 

iii) u satisfies for almost every x G R n 

iv ) J^T = A ^,g( u +) on fln > 
(±1 + S c A w , g - B c )f = e*- x on dQ 




(EPW 



(BE) /etf^an). 



Then all these problems are equivalent, in the sense that if a solution exists 
(is unique) for one problem, then a solution exists (is unique) for all the 
problems. If u is a solution of (DE), then u solves (IE), u\q> solves (EP), 
and u\qq solves (BE). 

Remark. The assumption that is not an eigenvalue is for simplicity. With 
appropriate changes, similar results are valid also when is an eigenvalue. 

The first step is to show the equivalence of the differential equation (DE) 
and the integral equation (IE). This involves the Green function G^, defined 
by 

G c = e**g c 

where g^ is the tempered distribution such that A^ 1 / = g^ * f for / in the 
Schwartz class (g^ exists since A^ 1 is translation invariant). Then 

AG C = C 2 G ? + 2e**C • Dg ( + e^ x Ag ( = A = S 



where So is the Dirac measure at 0. Consequently = Gq + where 

solution of 

n(n— 2)a(n) 



Gq(x) = c n \x\ 2 n is the usual fundamental solution of A, and Hr is a global 



harmonic function (one has c„ = —> — h ? \ where a(n) is the volume of the 

V " n(n— 2)a(n) V / 

n-dimensional unit ball). 

We note that the left hand side of (IE) is well defined for any u G 
Hl c (R n ), since then 2W ■ Du + {W 2 + D ■ W + q)u G L^(R n ) is a compactly 
supported distribution. Also note that G^ * f = e l< >' x A~^ 1 e~'' < '' x f whenever 

/ G L 2 M n )- 
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Lemma 5.1. Assume the conditions of Proposition 15.11 Then, u is a solu- 
tion of (DE) if and only if u is a solution of (IE). Also, a solution of (DE) is 
unique if and only if a solution of (IE) is unique. 

Proof. Let first u = e^' x (l + u>) solve (DE) where u> = A^ 1 / with / G Lq. 
Clearly u G Hy oc (R n ), and Hyy >q u = implies 

(A c + 2W ■ D c + (W 2 + D ■ W + q))(l + A" 1 /) = 0. 

We have A^(l + A^ 1 /) = /. Now applying A^ 1 to both sides, which is 
allowed since the left hand side is in Lq, gives 

u + A- 1 {2W • D f (1 + u) + (W 2 + D ■ W + q)(l + to)) = 0. 

We obtain (IE) by adding the constant one to both sides and multiplying 
by e^' x . 

For the converse, suppose u solves (IE), and write u = e^' x UQ. Then uq 
solves u + A^ 1 (2W-D (: u + (W 2 + D- W + q)u ) = 1. Applying A c to both 
sides gives B.-\^^ q u = 0. Also, one sees that u$ — 1 = A^ 1 / for / G L^(R n ). 

The uniqueness part is obtained just by noting that if u\ and solve 
(DE) then ui and U2 solve (IE), and vice versa. □ 

Next we show that (IE) and the exterior problem (EP) are equivalent. 
We have used the notation O' = R™ x f2 and Q' R = B(0,R) \ f2, where 
R > Rq and 0, C B(0,Rq). We write u+ (resp. for the restriction of u 
to 90 from the exterior (resp. interior), and (resp. g^-) for the value 
of Vu • v on <9$7 from the exterior (resp. interior), where v is the outer unit 
normal to dQ. We also write G^(x,y) = G^(x — y). 

A main point will be that a solution u of (IE) satisfies the radiation 
condition 

for a.e. x G R n as R — > oo. For applications of Green's identity below, we 
define a smooth approximation of G^ by G^ = Gq + where 

G e (x) = c n ( e 2 + |x| 2 )^. 
Note that AG^(x) = e~ n ip(x/e) where 

/ \ 1 /- i i?-. "+ 2 

^ =-rr 1+x 2 
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and J ip{x) dx = 1. Thus AG^ is an approximation of the identity. 

Before showing the equivalence of (IE) and (EP) we need a lemma on 
regularity properties of solutions of Hw, q u = and of Aw,q- 

Lemma 5.2. Under the conditions of Proposition 15.11 the operator Pw,q, 
which maps / E H 3 ^ 2 ({1) to the solution u of Hw q u = in £1 with u\qq = /, 
is bounded H 3 / 2 (dCl) -» H 2 (n). Further, one has A w>q : H 3 / 2 (dn) -» 
H^idn), and 

du 



&W,qf = q~ 

av 



an 



Proof. The operator Hw q , written in nondivergence form, satisfies the as- 
sumptions of Theorem 8.12] (the theorem is given for C 2 domains but 
the result holds with the same proof for C ♦ 1 domains). This shows that u is 
in H 2 (Q) if / E H 3 / 2 (d£l), and that the solution operator Pw,q is bounded. 

For the second part, we claim that if W E LQ(R n ;C n ) and D ■ W E 
L n/2( R n. c ^ then for any v G w i, n /{n-l)(ty one has 



f (W-Dv + (D- W)v) dx = 0. 
Jn 



(23) 



This statement means that W ■ v = on dQ, in a certain weak sense. The 
expression is well defined since v E L n /( n ~ 2 ) by Sobolev embedding. We 
take Wj E C£°(R, n ;C n ) to be convolution approximations of W so that 
Wj — > W in L n and D ■ Wj — > D ■ W in L n l 2 , and we take an extension of v 
in ^l.«/("-l)(R"). if the supports of and TU are contained in B(0,R), 
then 



/ (W-Dv + (D- W)v) dx = lim / {Wj ■ Dv + (D ■ Wj)v) dx 
Jn ^°°Jb(o,r) 

= lim - f (Wj ■ v)vdS = 0. 

i^oo l Jqb(o,R) 



Now let f,g£ H 3 / 2 (£l) and let Uf = Pw,qf and e 9 E H 2 (£l) with e 5 |aQ = 5. 
An integration by parts gives 

,<?)= /(Vn / -Ve 9 + (2W-L'u / + (W 2 + J D-W r + g)u / )e 9 )(ix. (24) 
Now u /e;? E TU 2 ' 1 ^) C T^Wfa-l)^). Using © with u = u/e 5 and 



substituting this to (|2*1)) gives 

is bounded H 3 / 2 (dn) -> H 1 / 2 ^). □ 



= Av^q/. This also shows that Ajy g 

90 
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Lemma 5.3. Assume the conditions of Proposition 15.11 Then, if u is a 
solution of (IE), then u\qi is a solution of (EP). Conversely, if u is a solution 
of (EP), then there is a unique extension u of u to R n so that u is a solution 
of (IE). Also, a solution of (IE) is unique if and only if a solution of (EP) is 
unique. 

Proof. Suppose u solves (IE). By Lemma 15.11 we have Hy/ )q u = and u G 
H 2 (H n ), which shows (EP) i)-ii). To prove hi) fix x G R n and let R > \x\ 
and R > Rq, and write 



/ du dGf(x,y)\ 

(Gl( X ,y)-(y)- U (y)-£gi)dS(y) 



\y\=R 

= - f (GUx, y)Au(y) - u{y)\G\{x, y)) dy 
Jb(o,r) 

= [ uA y G e Ax,y)dy+ [ G £ Ax,y){2W ■ Du + (W 2 + D -W + q)u) dy 

JB(0,R) JB(0,R) 

= (AG e ( * u X b(o,r))(x) + (Gl *(2W-Du + (W 2 + D-W + q)u))(x) 

since W and q have their supports inside -B(0, R). 

As e — * 0, the first term on the right converges to u{x) outside a set of 
measure zero (this set depends on R, but one may take the union of such 
sets for countably many R). The second term on the right converges to 
(G c * (2W ■ Du + (W 2 + D ■ W + q)u))(x) for a.e. x £ R ra by dominated 
convergence, since G^ £ L\ QC and the other function is in L 2 . Since \x\ < R 
the boundary integrals present no problem and one may replace G^ by G^. 
We obtain for a.e. x 



ihn / ( G({x , s A s) - uiy) d ^l) dSi9) 

R-.<x>./i„i = D \ dv dv(y) J 



R^ooJ lyl=R \ ' du dv(y) 

= u(x) + (G ? * (2W ■ Du + (W 2 + D • W + ?)w))(a;). (25) 

Since u satisfies (IE) we obtain (EP) hi). Further, since u G H? and is 
not a Dirichlet eigenvalue of Rw,q i n ^> Lemma 15.21 gives 

du du 



dv + dv- 
which is (EP) iv). 

Let now u solve (EP). We use Lemma and let v = P\v, q u+ G H 2 (Q), 
and we define u{x) = u{x) for x G Q' and n(x) = i>(x) for x £ Q. Now 
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u- = V- = u + = u + and 



du dv du du 



dv- dv- ' dv + dv + 

by Lemma ICT and (EP) iv). This shows that u G #£ c (R n ). By (EP) i) we 
have H\v )q u = in R n , and then the computation above leads to (|2*5|) with 
u replaced by u. The condition (EP) iii) shows that u solves (IE). 

The uniqueness part follows from the facts that if u\ and 112 solve (IE) 
then ui\qi and U2\w solve (EP), and if u\ and u<i solve (EP) then u\ and U2 
solve (IE). □ 

The final equivalence will be between (EP) and the boundary integral 
equation (BE). Here we need the layer potentials depending on f, defined in 
terms of the Green function G^. The single layer potential Sf, double layer 
potential D^, and boundary layer potential are defined by 

S ( f(x) = [ G c (x, y)f(y) dS(y) (x G R n \ M), 
Jan 

D ( f(x)= [ dG ^\ y) f(y)dS(y) (^R\ffl), 
Jan dv VV) 

Bd(x)= [ dG ^\ y) f(y)dS(y) (xedU). 
Jan 9u{y) 

Since d£l is C ' one does not need a principal value in the definition of B^. 
The operators have the following properties, given in |13j . 

Proposition 5.2. Let C R n , n > 3, be a bounded domain with C 1 ' 1 
boundary, and suppose that C B(0,Rq). 

(a) Let / G H l ' 2 {d£l) and u = Stf. Then u is in H 2 (Q) and H 2 (n' R ) for 
any R > i?o, and An = in R n \ d£l. If i? > Rq, then ii satisfies the 
radiation condition 

/^K'*^-^ 2 ^)^" (26) 

for almost every rc with |x| < it!. 

(b) Let / G H^l 2 (d£l) and u = D^f. Then v has the properties listed in 
(a). 



22 



(c) In the situation of (a), one has u_ = u+ on d£l, in the sense of 
H^l 2 (d£l) as well as nontangential convergence a.e. on <9f2. We will 
write u = S^f on d£l. The map / i— > S^f is bounded H s (d£l) — > 
H s+1 (dQ) for < s < 1, and one has 

^ ^ / an«l. (27) 



(d) In the situation of (b), one has 

v ± = ±^f + B c f ondn, (28) 

in the sense of H 3 / 2 (dQ) as well as nontangential convergence. 

(e) The map B c is bounded H s (dn) -» H s (dn) for < s < 3/2. 

Lemma 5.4. Assume the conditions of Proposition 15.11 Then, if u is a 
solution of (EP), then / = u\gn is a solution of (BE). Conversely, if / is a 
solution of (BE), then 



JC-x 



S ( A W J + D c f (29) 



is a solution of (EP), with u + = / '. Also, a solution of (EP) is unique if and 
only if a solution of (BE) is unique. 

Proof. Suppose u solves (EP). We let / = u + on dQ. Then / £ H 3 / 2 (dn). 
If x G and i? > 1x1, we have 



(G£(x, y)Au(y) - n(2/)A y G^(x, y)) cfy 

= (/ "/ ) (g £ c (^)|%) - u (y) dG } { (\ v) ) dS(y). 

\J\ v \=r Jdn JK ov dv(y) J 

Letting e — > and using (EP) i) we obtain 

" ( G <<^)f>>-«<> ^f 1 ) ^W-5 c (^)(x) +D{ („ + ) (l ) 

(30) 

for a.e. x in 17'. We let R — ► oo, use (EP) iii)-iv), and then let x — > 
<9f2 nontangentially and use Proposition 15.21 (d). which gives that / = u + 
satisfies (BE). 
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Conversely, suppose / satisfies (BE) and define u by (|2*5|) in Q! '. Then 
u satisfies (EP) i)-iii) by Proposition 15.21 (it is an easy calculation that the 
left hand side of $Zjj$ equals e l ^' x if u = e^' x ). We need to show (EP) iv). 
First note that by Proposition 15.21 (d). 

u+ = e^ x - S ( A W J + i/ + B ( f on Oft, 

which gives u+ = f using (BE). The formula (|30j) holds for u with the same 
proof, and as R — ► oo we obtain 

«(x) = e<--5 c (^-)(x)+2? f /(x) 

a.e. in f2'. Comparing with Q29JI we get 

Sc(|^-A^/)=0 (31) 

a.e. in O'. This holds also on dQ by Proposition 15.21 (c). and the uniqueness 
in the Dirichlet problem for A in fl shows that we have (|31|) in R n . Then 
(|27|) gives that jj± = A W J on dQ. 

If u\ and it2 solve (EP) then u\\qq, and «2|ar2 solve (BE), and if f\ and 
ji solve (BE) then the corresponding functions defined by (|29|) solve (EP). 
This shows the uniqueness part. □ 

Remark. Following Nachman [Tl], (BE) is equivalent to 

(/ + S c (A w>q - Ao l0 ))/ = e** on 0fi 
/ G H 3 / 2 {dQ). 

This follows since for x S fi' 

= (G c (x, .),Ao,o/) = 5 c (A ,o/)(x) 
and letting x — > c?J7 nontangentially gives ^1 + = S^Aq^. 

Proposition 15.11 is an immediate consequence of Lemmas 15.11 to 15.41 We 

conclude the section by showing that the operator arising in (BE) is of 
the form I + K with K compact. This fact and the Fredholm alternative 
show, for instance, that uniqueness in one of the problems in Proposition 
15.11 implies the existence of a unique solution for all the problems. 
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Lemma 5.5. Let f2 C R n , n > 3, be a bounded domain with C ' boundary. 
Then the operator S c A W:q - B c - \l : H 3 / 2 (dVt) -> H 3 / 2 (dn) is compact. 

Proof. Let / 6 # 3/2 (dft) and let u = If a? € fi we have 

- / (^(x,y)Au(y)-u(y)A y G^,y))(iy 
./n 



If e — > we get 

u(ar) + / G c (x, y)(2W ■ Du + (W 2 + D ■ W + q)u) dy = (S ( A w , q - D ( )f(x) 
Jn 

a.e. in f2. Let then x — > dQ nontangentially, so that Proposition 15.21 gives 

{S c K w , q -B c -\l)f = R [ G ( (x,y)(2W-D + (W 2 + D-W + q))P w , g f(y)dy. 
1 Jn 

This reads 

S ( A W)g -B c -h = RG c MJP w>q 

where R is the trace H 2 {ti) -> H 3 / 2 (dn), G ( : L 2 (Q) -> F 2 (ft) restricts 
G^*u = e 1 ^'* A.^ 1 e~^' x u to where it is the extension by zero of u £ L 2 (Q) 
to R n , M : H 1 ^) -> L 2 (ft) maps u to 2W • + (W 2 + L> • W + g)u, and 
J is the embedding F 2 (S1) -> iT 1 ^). All these maps are bounded and J is 
compact, so the composition is compact. □ 



6 Reconstruction of the magnetic field 

The preceding section considered equivalent formulations for problems which 
give rise to CGO solutions, but did not consider the solvability of any of the 
problems. The next result, which follows directly from Theorem ll.il shows 
that if W is continuous then the problems indeed have unique solutions for 
large £. 

Proposition 6.1. Assume the conditions in the beginning of Theorem EH 
Then there exists C = C(n, f2, W, q) so that whenever \(\ > C, then each of 
the problems (DE), (IE), (EP), (BE) has a unique solution. 
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Proof. It is enough to show that (DE) has a unique solution. Now u = 
e l< *' x (l + uj) solves Hy/ A u = if and only if 

(A ? + 2W • D ( + (W 2 + D ■ W + q))u = -(2£ • W + W 2 + D • W + q). (32) 

From Theorem 11,11 we know that if |£| > C(n, J7, W, q) this equation has a 
unique solution uj G i^j . Then l|3*2^1 gives that A^uj = f for some / G Lq , so 
that cj € A^L^. This shows that u is the unique solution of (DE). □ 

We will from now on assume the conditions in the beginning of Theorem 
11.21 For given £ we denote by the unique solution of (DE). It follows that 
if one knows h-w,q then the boundary values u^\qq may be reconstructed as 
the unique solution of the boundary integral equation (BE) . The rest of the 
section will be devoted to showing that the magnetic field curlPF may be 
reconstructed from this knowledge. The first step, similarly as in [13., is to 
consider a (non-physical) scattering transform. 

Definition. Let ^ G R n be such that |£| 2 is not a Dirichlet eigenvalue of A 
in 0. Then for any ( G C n which satisfies C 2 = 0, |C| > C, ReC -L f, and 
Im £ _L £, we define 

It is clear from the preceding discussion that h.w,q determines tw,q for 
the appropriate £, £. Using the weak formulation of the DN map and the 
fact that i?o,-|£|2e - * x '^ + ^ = in $7, one easily sees that 

tw,q&Q= [ e- ix <(2((-W)uo+W-Du + (W 2 H-W+\Z\ 2 +q)u )dx, (33) 
Jn 

where we write uq = e~ l< ^' x u^. 

We know from (DE) that no = 1 + uj with uj G A^ Lq. If one had 
IMIl 2 (q) — ► as |£| — > oo one could divide (|3"3*|) by \(\ and let \(\ — > oo, 
which would then give essentially the Fourier transform of curlPF. However, 
uj is obtained by solving (|32|) where the L 2 +1 norm of the right hand side is 
O (|C|) instead of o(|£|), so one gets that ||<-^||z, 2 (ri) is bounded but may not 
be small when |£| is large (in fact Lemma Id II shows that uj — * e l< ^ — 1 in 
L 2 (Q) as \C\ — ► oo, where 4> is defined below). 

To deal with this difficulty we write the solution in a different form 
where one gets a small remainder term for |£| large. For this we employ a 
decomposition 

W = W i + 
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where = W * cp e with ip e the usual mollifier, and we make the choice 

e = icr 

with < a < 1/2. Then is C°° and 

\\W*\\ w i^=o{\C,n (34) 
111^11^-= (|C| 2 '), (35) 
||^|| L oo = o(l) (36) 

as |C| — > oo. 

We also write £ = sfi where /U = 71 + 272, |7j| = 1, 71 -1-72- Finally, we 
fix X e C c °°(R n ) with x = 1 in B(0,M/2), x = outside of B(0,M), and 
H C £(0,M/2). 

Lemma 6.1. Fix > with a + < 1/2. For \(\ large enough, the CGO 
solution of Hwqii^ = in R" may be written in the form 

u c = e^ x {uj + uj) (37) 

where ojq = e ix ^ with 

X( ( x )= x ( x /\Cf), (38) 

0»(x)=iV / 7 1 (- / u.^), (39) 

and to G ili with [|w||ra = o(l), ||u;||#i = o(|£|) as |£| — > 00. 

Proof. We first show that the equation Hw,qU = in R ra has a solution of 
the form (|37|) with the required properties. This will be the case if ui satisfies 

(A c + 2W • D ( + G> = -/ (40) 
where G = W 2 + D ■ W + q € and 



/ = (A c + 2VF • D c + G)u = e lx ^ ixc A ^ + 2iD X( ■ Dfl + i^A X<: 
+ (x C V0« + ^Vxc) 2 + 2C • (V XC )^ + 2C • (V0«)xc 

+ 2W • (Vxc)0 tt + 2VF • (V0 tt )xc + 2VF fl • C + 2W b • C + G 

We need to know the behaviour of ||/||r,a as |CI grows. The choice of 

(5+1 

implies 

2C • V(f)U + 2W 9 ■ C = 0. 
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bmce = X(Wt this removes the worst two terms from /, and one obtains 



in terms of L| . j norms 



< C 



HxcA^II + livxc • V0»|| + II^Axdl + lllxcV^I 2 !! + |||^vxcl 2 ll 
+ |cMl(Vx(*/|Cl e ))^ll + \\w ■ (Vxd^W + \\w ■ (v^)xdl 

+ ||^. C || + || G ||]. (41) 

Lemma 13.11 implies 

\d a 4>Kx)\ < cicr^wrWi^) (42) 

where xt is the projection of x to span{7i, 72} and x± = x — xt- Then for 
instance 

llxcA^H = (j^(x) 2 ^ X<: (x) 2 \A^(x)\ 2 dxy /2 

<C\(\ 2 °( f (x) 2 ^(x T )- 2 dxY /2 

^ J\x T \<M\(\ e ,\x ± \<M ' 

< C\C\ 2a ( [ (x T ) 2S dx T Y /2 < C \{\ 2 °+( s +V e . 

^ J\x T \<M\C\ e ' 

This has the worst behaviour of the first five terms of (|41|) since derivatives 
hitting x<; bring decay in |£| and the other terms involve only first derivatives 
of A similar computation shows that the sixth term is 0(|£| 1 ~ e +( <5 + 1 ) 6 ') = 
0(\(\ 1+Sd ). One has W ■ Vxc = and Wx( = W for large |C| so the seventh 
and eight terms are and 0(|£| CT ), respectively. The final two terms are 
o(|C|) and O(l), respectively, since \\W as |C| — ► 00. Using the 

choices of a and 6 and the fact that — 1 < 5 < 0, we obtain ||/|| L i = o{\Q\). 
The solution to of (|40jl has the desired properties by Theorem ll.il 

It remains to show that u given by (|37j) is the CGO solution. One has 
u = e l ^' x (l + Cj) where 

Q = jX^ _ I + w> 

Now e lx ^ — 1 is in Hj for any 5 < 0, since for instance e lx ^ — 1 = 0(\x( 
and (|42|) implies 



||e*c*» - < C(^(x) 2<5 xc (x) 2 |^(x)| 2 dx) 1/2 
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Also lo G Hi so uj G i/r for — 1 < 5 < 0. This and Hy/^ q u = imply 
u) G A^Lq, so u is indeed the unique solution of (DE) given by Proposition 
EH □ 

We may now plug in uo = ujq + u> from lpT7|) to . The estimates for 
u and the form of ujq imply that 

= lim s-^woC^sm) = 2 / e -<x V*(p ■ W) dx 

s^co " J 

where 4> = N~ 1 {—^l ■ W). This shows that we may recover a nonlinear 
Fourier transform Rw,q(^, A*) of // • W from the knowledge of Aw,q, for any 
[i and £ with £ • \i = and |£| 2 not a Dirichlet eigenvalue of A in f2. 

The next argument, due to Eskin and Ralston 5 , shows that this nonlin- 
ear Fourier transform is in fact just an ordinary Fourier transform. Similar 
ideas appear in Sun [2*3] . 



Lemma 6.2. One has 



Rw, q (Z> A*) = 2 y e- iX - 5 (//-Ty)(ix. 



Proof. It is enough to prove this for \i = e\ + ie 2 , so that £ = (0, £") and 
(9i + id 2 )4> = -(Wi + iW 2 ). Then 



R W ^,V) = 2 / e"" * e i ^{-{d l +id 2 )ct>)dx 

JR" 

= 2/ e-^'^'/i^'Odx" 



R 



n — 2 



where 



/i(x") = * / (fli + id 2 )(e i ^ x > x ' } ) dx' = lim t / (Si + ^)(e^ ')) dx' 

iR 2 fl - >0 ° J\x'\<R 

lim t / e^'^'Vi + ^^x'). 



R—*oo 



x'\=R 



The integrals are well defined by standard approximation arguments. Now 
e *0 = 1 + i<j, + <3(|#| 2 ) = 1 + i<j> + 0(|x'|- 2 ) by Lemma ETT1 and 

J lx , l=R (vi + ^2) dS(x') = f lx ,^ R (di + id 2 )(l) dx' = 0, 



L,i = „o(ixr 2 )(^i+^2)^(x') 



< § — > as — > 00, 
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so we have 



R— >oo 



h{x") = - lim / 4>{x' ,x"){vx + iv 2 ) dS{x') 



x'\=R 



= — lim / (di + id2)4>dx' = (e\ + ze2) • W dx' . 

R-^OO J\ X >\< R J R 2 

This gives the claim. □ 

We now show that Rw,q{£,,lj) determines curlVF, or DjWk — DhWj for 
any j ^ k. Let £ G R n be such that |£| 2 is not a Dirichlet eigenvalue of A 
in f2. If one of ^ is nonzero choose 7 to be the unit vector with direction 
^e/c — 6c e i) so that £ • 7 = 0. Since n > 3 we may choose a unit vector 
7 with 7 • 7 = £ • 7 = 0. Letting /j = 7 + 17 we see that Aw^ determines 
^W,?(£>/f) + -Rw,g(^)A*)) which determines 

J e- lx <(CjW k - tkWj) dx = (DjWk - DkWjYiO 

in fact for any £ G R n such that |£| 2 is not an eigenvalue. Since there are 
countably many eigenvalues and since DjWk—D^Wj is compactly supported, 
so the Fourier transform is analytic, we recover DjWk — D^W j. 



7 Reconstruction of the electric potential 

Finally, we make the additional assumptions that W is C 1+e and d£l is 
C 2+£ for some e > 0, and we indicate how to recover q from A.w t q- From the 
preceding section, we may assume that the magnetic field curl W is known. 
The next step is to construct a certain magnetic potential with this magnetic 
field. 

Lemma 7.1. One can construct W G C 1+e (f2;C n ) with curll^ = curlVF 
and W\ an = 0. 

Proof. For the following concepts we refer to [U- We write XA k (Q) for 
a differential A;- form in Q with coefficient functions in X, and tr\ and rvq 
for the tangential and normal traces on <9f2, respectively, of a form 77. Let 
ijj = Y^j=i Wjdxj £ C 1+£ A 1 (Q) be the 1-form corresponding to W, and let 
X = duj. We start with solving the boundary value problem 

f dCb =x in 0, U o\ 
}&=0 onffi { ' 
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for u) in C 1+£ A 1 (Q). In the case of smooth domains and L p Sobolev spaces 
this problem is considered in j^J , where the solution is reduced to the Hodge 
decomposition of \. For a general form \ S L 2 A k (Q) this decomposition 
reads 

x = da + 5(3 + k (44) 

where a £ ffU^^) with to = 0, (3 £ ff^^^fi) with n(3 = 0, and k 
is a harmonic field meaning that dn = 5k = 0. Here 5 is the codifferential. 
Further, the three summands in (|44j) are uniquely determined and mutually 
orthogonal with respect to the natural L 2 inner product. 

The specific form of \ above implies that the decomposition ()44l) reduces 
to x = da, where one may choose a = 5<p where <j) is the Dirichlet potential 
of x (see Section 2.2). Under the present assumptions of C 2+e boundary 
and C £ regularity of x> Theorem 7.7.4 in |12j implies that 5(p £ C 1+£ . It 
is then easy to check that lv = 5<j) solves (|43[) , We note that the Dirichlet 
potential <fi is obtained constructively using the explicit integral formula for 
the corresponding Green operator, as in 

Letting W be the vector field corresponding to a), we have that W £ 
(7 1+e (0; C ra ), curlPU = curl W, and the tangential components of W vanish 
on dCl. We may further replace W by W + Vp where p £ C 2+e {Q) satisfies 
p\d£i = and ^jldfj = —W-u, and p is constructed similarly as in Lemma 
5.8 of [201 • This completes the proof. □ 

With W as above, we conclude that W = W + Vp where p £ H /1,00 (0) 
and p\on = 0. An easy argument for this under the present regularity 
assumptions is obtained by extending W — W by zero to R™ as a Lipschitz 
vector field, and by noting that W — W = Vp in a large ball with p £ W 2 '°° 
and Vp\gn = 0. The assumption on the topology of 17 ensures that 8ft is 
connected, so that p is constant on dfl, and one may substract the constant 
to get p\on = 0. Gauge equivalence then implies that = Aw,o- 

Fix £ £ R n n {0} and take unit vectors 7,- such that {£,71,72} form an 
orthogonal set. For s > define complex vectors 

Cl = _ 2 +S V 4^ 71+ ' S72 ' 

(2 = ~2~ s y l ~A? ii ~ tsi2 - 

Using the notation of Lemma KTT1 the equations Hw, q u = and H_w$v = 
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have unique CGO solutions u = and v = v^ 2 , which have the form 

t, Ca = e * a -*( e -*c a ^+w 2 ) 

where 11^,2 — > as s — > oo. 

We define a new scattering transform 

*(0 = ((Aw, g - a^o)^!^),^!^). 

Since A^q an d A_^q are known, one may construct the boundary values 
of 11^ and v^ 2 as solutions of boundary integral equations as in Section 5, 
and thus t is known. The definition of DN maps implies 

t(0 = J Q u (i v (2 dx, 

and so — > J e^^q^x) as s — > oo. This is the Fourier transform of q, 
and we have recovered the electric potential. 
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